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We derive a formula for the entropy for a multicomponent coupled fluid, which under special 
conditions reduces to the Cardy-Verlinde form relating the entropy of a closed FRW universe to its 
energy together with its Casimir energy. The generalized fluid obeys an inhomogeneous equation 
of state. A viscous dark fluid is included, and also modified gravity is included in terms of its fluid 
representation. It is demonstrated how such an expression reduces to the standard Cardy-Verlinde 
formula corresponding to the 2d CFT entropy in some special cases. The dynamical entropy bound 
for a closed FRW universe with dark components is obtained. The universality of the dynamical 
r-\ • entropy bound near a future singularity (of all known four types), as well as near the Big Bang 

_^r* singularity, is investigated. It is demonstrated that except from some special cases of Type II and 

Type IV singularities the dynamical entropy bound is violated near the singularity even if quantum 
C^*) ■ effects are taken into account. The dynamical entropy bound seems to be universal for the case of 

CNJ ' a regular universe, including the asymptotic de Sitter universe. 
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I. INTRODUCTION 



Recent astronomical data indicate that the current universe is expanding with cosmic acceleration caused by the 
so-called Dark Energy (DE). The well-known ACDM model where Dark Energy comes from the effective cosmological 
constant fits quite well the observational data. Nevertheless, it is not excluded that the current DE may have a 
— — ■ phantom origin (or will become phantom- like DE in the near future). Moreover, the possibility of a quintessence- like 
DE with effective equation of state (EoS) parameter being very close to —1 is widely discussed in the recent literature. 
It is quite likely that the current FRW universe is spatially-flat. However, the occurrence of a slightly spatially-curved 
' FRW universe is not excluded by observational bounds. 
t-H . The study of phantom-like or quintessence-like effective dark fluids opens for a number of new phenomena which 
are typical for such a DE universe. For instance, it is known that phantom DE may drive the future universe to a 
. rH ! so-called finite-time Big Rip singularity (for earlier works on this, see Refs.[l], [2[). From another side, quintessence-like 
DE may bring the future universe to a milder future singularity (like the sudden singularity ,3;, 4j where the effective 
energy-density is finite) . Actually, the study of Ref. [|[ shows that there are four different types of future finite-time 
singularities where the Type I singularity corresponds to the Big Rip, the sudden singularity is of Type II, etc. The 
universe looks quite strange near to the singularity where curvature may grow up so that quantum gravity effects 
may be dominant In any case, the study of the universe under critical conditions (for instance, near a future 
singularity) may clarify the number of fundamental issues relating seemingly different physical theories. 

Some time ago [f| it was shown that the first FRW equation for a closed FRW universe may have a more fundamental 
origin than what is expected from standard General Relativity. It was demonstrated that this equation may be 
rewritten so as to describe the universe entropy in terms of total energy and Casimir energy (the so-called Cardy- 
Verlinde (CV) formula). Moreover, it turns out that the corresponding formula has a striking correspondence with 
the Cardy formula for the entropy of a two-dimensional conformal field theory (2d CFT). Finally, the formula may 
be rewritten as a dynamical entropy bound from which a number of entropy bounds, proposed earlier, follow. The 
connection between the standard gravitational equation and the 2d CFT dynamical entropy bound indicates a very 
deep relation between gravity and thermodynamics. It raises the question about to which extent the CV formula is 
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universal. Problems of this sort are natural to study when the universe is under critical conditions, such as near a 
singularity. 

The present work is devoted to a study of the universality of the CV formula and the corresponding dynamical 
entropy bound in a DE universe filled with a generalized fluid, especially near the singularity regime. Generalization 
of the CV formula for a multicomponent fluid with interactions, assuming the EoS to be inhomogeneous, is presented. 
The viscous case is incorporated. It is shown that the standard CV formula with correct power (square root) is 
restored only for some very special cases. The dynamical entropy bound for such fluids near all the four types of 
the future singularity is considered. It is demonstrated that this dynamical entropy bound is most likely violated 
near the singularity (except from some cases of Type II and Type IV singularity). This situation is not qualitatively 
changed even if account is taken of quantum effects in conformally invariant theory. Using the formalism of modified 
gravity to describe an effective dark fluid, the corresponding CV formula is constructed also for i 7 '(i?)-gravity. The 
corresponding dynamical entropy bound is derived. It is shown that the bound is satisfied for a de Sitter universe 
solution. Further discussion and outlook is given in the discussion section. 



II. GENERALIZATION OF CARDY-VERLINDE FORMULA IN FRW UNIVERSE FOR VARIOUS 

TYPES OF FLUIDS 

This section is devoted to consideration of the Cardy-Verlinde formula for more general scenarios than those 
considered in previous works (see @-[3|)- We consider a (n + l)-dimensional spacetime described by the FRW metric, 
written in comoving coordinates as 

where k = — 1, 0, +1 for an open, flat, or closed spatial Universe, a(t) is taken to have unit of length, and dil^-i is the 
metric of an n — 1 sphere. By inserting the metric (JT]) in the Einstein field equations the FRW equations are derived, 



o 16irG v-^ k ■ 8nG v-^. , k 

H 2 = - ttI>--' H = j"I>*+J»i) + -2 ■ 2 

i—l 



n(n - 1) ^ 



Here pi = Ei/V and pi are the energy-density and pressure of the matter component i that fills the Universe. In 
this paper we consider only the k = 1 closed Universe. Moreover, we assume an equation of state (EoS) of the form 
Pi — WiPi with Wi constant for each fluid, and assume at first no interaction between the different components. Then, 
the conservation law for energy has the form 

p\ +nH(pi +pi) = , (3) 

and by solving ([3]), we find that the i fluid depends on the scale factor as 

Pl cx a ^ 1+w ^ . (4) 

Let us now review the case of Ref.Q, where just one fluid with EoS p = wp and w = constant is considered. The total 
energy inside the comoving volume V , E = pV, can be written as the sum of an extensive part E E and a subextensive 
part Eq, called the Casimir energy, and takes the form: 

E{S,V)=Ee(S,V) + ±E c (S,V) . (5) 

Under a rescaling of the entropy (S — > AS 1 ) and the volume (V — > XV), the extensive and subextensive parts of the 
total energy transform as 

E E (XS,XV) = XE E {S,V) , E C {XS,XV) = X^^EdS, V) . (6) 

Hence, by assuming that the Universe satisfies the first law of thermodynamics, the term corresponding to the Casimir 
energy Ec can be seen as a violation of the Euler identity according to the definition in Rcf . [6] : 

E c =n(E + pV-TS) . (7) 

Since the total energy behaves as E ~ a~ nw and by the definition the Casimir energy also goes as Eq ~ a~ nw . 
The FRW Universe expands adiabatically, dS — 0, so the products Eca nw and E E a nw should be independent of the 
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volume V, and be just a function of the entropy. Then, by the rescaling properties (J6j) , the extensive and subextensive 
part of the total energy can be written as functions of the entropy only 0] , 



E E = -^S w+1 , E c = ^—s w+1 ' 2 / n . (8) 

Here a and j3 are undetermined constants. By combining these expressions with ([5]), the entropy of the Universe is 
written as a function of the total energy E and the Casimir energy Ec , [7] , 
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S = ^-^=-^E c (2E-E c )j , (9) 

which for w = 1/ n (radiation-like fluid) reduces to @ , 

S=^=VE C (2E-E C ) , (10) 

which has the same form as the Cardy formula given in Ref.[l]. The first FRW equation @ can be rewritten as a 
relation between thermodynamics variables, and yields 

Sh = —a^E BH {2E-E BH ) , where S H = (n-l)— , E BH = n(n - 1)-—^ . (11) 

It is easy to check that for the bound proposed in Ref.[(|, Ec < E BB , the equation for the entropy (fTUl) coincides 
with the first FRW equation (fTTj) when the bound is reached. We will see below that when there are several fluid 
components, the same kind of expression as in Ref.[6j cannot be found. Nor is there the same correspondence with 
the FRW equation when the bound is saturated. 



A. Multicomponent Universe 

If m fluids are considered with arbitrary EoS, pi — WiPi, the expression for the total entropy is simple to derive 
just by following the same method as above. The total entropy is given by the sum of the entropies for each fluid, 



/ n wn nwi \ „ (tl , i + i)_i 

5 = E s * = E \—m~^ EiC ^ Ei - E ^)) ■ (12) 

This expression cannot be reduced to one depending only on the total energy unless very special conditions on the 
nature of the fluids are assumed. Let us for simplicity assume that there are only two fluids with EoS given by 
Pi — wipi and P2 — W2P2, Wi and W2 being constants. We can substitute the fluids by an effective fluid described by 
the EoS 

, P1+P2 w 2 -wi 

Pea = w cS p c fi , where w eS = ■ = W\ + — -. — , (13) 

P1+P2 I+P1IP2 

and p c ff = \{p\ + P2), Pes = \{p\ + P2)- Then, by using the energy conservation equation (|3]), we find p\ ~ 
{a/a )- n(1+wi) and p 2 - {a/a )- n(1+W2 \ where ao is assumed to be the value of the scale factor at the time t$. The 
effective EoS parameter w e g can be expressed as a function of the scale factor a(t) 

W2 — W1 . . 

W e B = Wi H ; — ; —, r . (14) 

The total energy inside a volume V becomes 

E T = E 1 + E 2 oc {a/a y nwi + {a/a Q )- nW2 . (15) 

As the energy is proportional to two different powers of the scale factor a, it is not possible to write it as a function 
of the total entropy only. As a special case, if the EoS parameters are W\ = w% = w c g, the formula for the entropy 
reduces to ([9]), and coincides with the CV formula when w e g = 1/n. 
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As another case, one might consider that for some epoch of the cosmic history, w\ ^> W2- Taking also a >> ao, 
we could then approximate the total energy by the function Ej< oc a~ nW2 . From ([5]) the Casimir energy would also 
depend on the same power of a, Eq oc a~ nW2 . The expression (j9|) is again recovered with w = w 2 - 

Thus in general, when a multicomponent FRW Universe is assumed, the formula for the total entropy does not 
resemble the Cardy formula, nor does it correspond to the FRW equation when the Casimir bound is reached. It 
becomes possible to reconstruct the formula (|10[) . and establish the correspondence with the Cardy formula, only if 
we make specific choices for the EoS of the fluids. 



B. Interacting fluids 

As a second case we now consider a Universe, described by the metric ([T]), filled with two interacting fluids. One 
can write the energy conservation equation for each fluid as 

pi + nH(pi + pi) = Q , p2 + nH(p 2 + p 2 ) = -Q , (16) 

where Q is a function that accounts for the energy exchange between the fluids. This kind of interaction has been 
discussed previously in studies of dark energy and dark matter. The effective EoS parameter is given by the same 
expression (IT51) as before. With a specific choice for the coupling function Q, the equations P^|) may be solved. One 
can in principle find the dependence of the energy densities p\^ 2 on the scale factor a, 

pi = a{t)- n{1+W1 ^ (d + J a n{1+Wl ^Q(t)d?j , p 2 = a{t)~ n< - 1+w ^ ( C 2~J a n ^ 1+w ^ Q{t)dt\ , (17) 

where C\ and C 2 are integration constants. In general it is not possible to reproduce the CV formula, and the result 
will be a sum of different contributions, similar to the entropy expression given in (|12[) . However, for the case where 
the effective EoS parameter (fT4"]) is a constant, the expression for the entropy will be given by equation © as before. 
This condition only holds when w c s = w\ = w 2 , where the situation is thus equivalent to the one-fluid case, and the 
entropy reduces to the CV formula when w c g — 1/n. 

Let us consider a simple choice for the function Q that leads to the CV formula for a certain limit. Let Q = Qoa m H, 
where m is a positive number, Q Q is a constant, and H(t) the Hubble parameter. Then the integral in (1171) is easily 
calculated, and the energy densities depend on the scale factor according to 

p x = C ia - n{1+w ^ + ha m , p 2 = C 2 a- n{1+w ^ + k 2 a m , (18) 

where kip — Qo/(n(l + W1.2) + to). If we restrict ourselves to the regime where a 3> Ci,2 such that the first terms in 
the expressions for p\^ 2 are negligible, the effective EoS parameter becomes 

w eS = wi + — T— . (19) 

Then, the entropy of the universe is given by ([9]) with w = w G ff . The CV formula can be reproduced only with very 
specific choice of the free parameters, just as above. 

We have thus shown that in general a formula for the entropy of the type ((9]) cannot be reconstructed for interacting 
fluids. Coincidence with the Cardy formula is obtained if the effective EoS parameter is radiation-like, u> c ff = 1/n. 
Then the expression for the entropy turns out to be in agreement with the formula (|10[) . corresponding to the first 
FRW equation (fTTj) when the Casimir energy reaches the bound Ec = £"bh- 



C. Inhomogeneous EoS fluid and bulk viscosity 

Let us now explore the case of an n + 1-dimensional Universe filled with a fluid satisfying an inhomogeneous EoS. 
This kind of EoS, generalizing the perfect fluid model, has been considered in several papers as a way to describe 
effectively the dark energy (see @ and [13]). We assume an EoS expressed as a function of the scale factor, 



p = w(a)p + g(a) 



(20) 
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This EoS fluid could be taken to correspond to modified gravity, or to bulk viscosity (Ref.|j[). By introducing (l20l) in 
the energy conservation equation (J3|) wc obtain 

n(l + w(a)) s(a) 

p (a) H p(a) = —n . (21) 

a a 

Here we have performed a variable change t = t(a) such that the prime over p denotes derivative with respect to the 
scale factor a. The general solution of this equation is 



p(a) = e- F ^ U-nJ e F ^ ^da^j where F(a) = f 1 ± J K) da' , 



(22) 



and K is an integration constant. As shown above, only for some special choices of the functions w(a) and g(a), the 
formula (|10[) can be recovered. Let us assume, as an example, that w(a) = —1 and g(a) — — a m , with m = constant. 
Then, the energy density behaves as p oc a m . Hence, by following the same steps as described above, the extensive and 
subextensive energy go as a m+n , and by imposing conformal invariance and the rescaling properties ([5]), we calculate 
the dependence on the entropy to be 

Eb = ~ A %—r -,S- m l n , E c = ^^o-(2n+W«) . (23) 

The expression for the entropy is easily constructed by combining these two expressions and substituting the 
extensive part by the total energy. This gives us the same expression as in ([9]) with w = — (n + m)/n. Note that 
for m — — (1 + n), the formula (|10[) is recovered and also its correspondence with the CFT formula. However for 
a generic power m, the CV formula cannot be reconstructed, like the cases studied above. Only for some special 
choices does the correspondence work, leading to the identification between the FRW equation and the Cardy formula. 

Let us now consider an inhomogeneous EoS fluid due to bulk viscosity. From a hydrodynamical perspective it is 
natural to extend the formalism so as to incorporate viscosity effects. Working to the first order in the deviations from 
thermal equilibrium we are faced with two viscosity coefficients, namely the shear viscosity r\ and the bulk viscosity £■ 
In conformity with spatial isotropy we shall assume, as usual, that only the bulk viscosity contributes. (For a review 
of viscous cosmology and entropy, one may consult Ref.[ll| and also Refs.[l2j]-) The viscous fluid may be considered 
as a special kind of inhomogeneous EoS fluid, although of a different kind from that of Eq. (|20|) above. We set the 
number of spatial dimensions n equal to 3. The energy- momentum tensor can be written as 

= pU^U v + ph^ , (24) 
where /i M „ = g M „ + U^U^ is the projection tensor and 

p = p-3H(, (25) 

the effective pressure. For simplicity let us assume in this subsection the following simple fluid model: 

w = constant, 5 = 0, £ = constant, (26) 

In comoving coordinates, U° = l,U l = 0. The first of the FRW equations ([2]) maintains its form (it is viscous 
insensitive), whereas the second equation becomes 

H = -47rG(p + p) + ^ , (27) 

showing that p is now the thermodynamically important pressure. We see that the relationship p — wp, or 

p = wp - 3H( , (28) 

can be considered as an EoS in the present case. 

Introduction of a viscosity means effectively the introduction of a length parameter, and so the conformal invariance 
of the formalism is lost. The question arises: Can the entropy arguments leading to the Cardy- Verlinde formula be 
carried over to the viscous case? The answer actually turns out to be affirmative, at least when £ is small. The most 
delicate point in the line of arguments is the assumed pure entropy dependence of the product Ea 3w . Now, if one 
uses the FRW equations to derive the "energy equation" (k = +1 assumed) 



di 



(pa z{1+w ^ = 9CH 2 a 3(1+tu) , (29) 
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one can combine this with the equation for entropy production 

Na = -±H 2 , (30) 

where N is the particle (baryon) density and a the entropy per particle. As discussed in some detail in [llj it follows 
that, since the total energy E ~ pa 3 and the total entropy S ~ Ncra 3 , the quantity Ea 3w becomes independent of 
the volume V and is a function of S only. This generalizes the pure entropy dependence of the product Ea found 
by Verlinde in the case of a non-viscous radiation dominated universe. We obtain the expression ((9J with n = 3 as 
the generalized Cardy- Verlinde formula and reducing to the standard formula (with square root) when the universe 
is radiation dominated. 

The following point ought finally to be noted. The energy conservation equation T 0l/ ; „ = implies 

p + 3H(p + p) = 0, (31) 
so that equations (21) can be taken over to the viscous case, only with the substitutions pi — >• pi. 



III. ON THE COSMOLOGICAL BOUNDS NEAR FUTURE SINGULARITIES 

InRef.0, Verlinde proposed a new universal bound on cosmology based on a restriction of the Casimir energy Ec] 
cf. his entropy formula (|10j) . This new bound postulated was 

E c < E Bli , (32) 

where -Ebh = n(n — 1) g^Ga 3 • ^ was deduced by the fact that in the limit when the Universe passes between 
strongly and weakly self-gravitating regimes, the Bekenstein entropy 5b = ^j^E and the Bekenstein-Hawking entropy 
<Sbh = (n — 1)4^1 which define each regime, are equal. This bound could be interpreted to mean that the Casimir 
energy never becomes able to reach sufficient energy, Ebh > to form a black hole of the size of the Universe. It is easy 
to verify that the strong (Ha > 1) and weak (Ha < 1) self-gravity regimes have the following restrictions on the total 
energy, 

E < £ B h for Ha < 1 , 

E > £ B h for Ha > 1 . (33) 

From here it is easy to calculate the bounds on the entropy of the Universe in the case when the Verlinde formula 
(fTUjl is valid; this is (as shown in the above sections) for an effective radiation dominated Universe w c g ~ 1/n. The 
bounds for the entropy deduced in Ref.Q for k = 1 are 

S < S B for Ha < 1 , 

S > S H for Ha > 1 , (34) 

where Sb is the Bekenstein entropy defined above, and Sn is the Hubble entropy given by (fTT|) . Note that for the 
strong self-gravity regime, Ha > 1, the energy range is Ec < £^bh < E. According to the formula (|10|) the maximum 
entropy is reached when the bound is saturated, Eq = E-qu- Then S = Sn, such that the FRW equation coincides 
with the CV formula, thus indicating a connection with CFT. For the weak regime, Ha < 1, the range of energies 
goes as Eq < E < Ebh and the maximum entropy is reached earlier, when Ec = E, yielding the result S = Sb- 
The entropy bounds can be extended to more general cases, corresponding to an arbitrary EoS parameter w. By 
taking the bound (I32p to be universally valid one can easily deduce the new entropy bounds for each regime, from 
the expression of the entropy (j9)). These new bounds, discussed in Ref.[7|, differ from the ones given in (|34"]) . but still 
establish a bound on the entropy as long as the bound on E c expressed in (|3"2"j) is taken to be valid. The entropy 
bounds can be related through the first FRW equation, yielding the following quadratic expression (for k = 1), 

<Sh + ($b — <Sbh) 2 = S^ ■ (35) 

We would like to study what happens to the bounds, particularly to the fundamental bound (|32[) . when the cosmic 
evolution is close to a future singularity; then the effective fluid dominating the cosmic evolution could have an unusual 
EoS. As shown below, for some class of future singularities such a bound could soften the singularities in order to 
avoid violation of the universal bound (j32[) . It could be interpreted to mean that quantum effects become important 
when the bound is reached. However, as the violation of the bound could happen long before the singularity even in 
the presence of quantum effects, it could be a signal of breaking of the universality of the bound (13"2"|) . Let us first of 
all give a list of the possible future cosmic singularities, which can be classified according to Ref.[5] as 
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• Type I ("Big Rip"): For t — > t s , a — > oo and p —> oo, \p\ —> oo. 

• Type II ("Sudden"): For t t s , a a s and p p s , |p| -> oo. (see Refs.0, 01) 

• Type III: For t — > t s , a — ¥ a s and p —¥ oo, \p\ — > oo. 

• Type IV: For t —>t s , a — > a s and p — > /9 S , p — > p s but higher derivatives of Hubble parameter diverge. 

Note that the above list was suggested in the case of a flat FRW Universe. As we consider in this paper a closed 
Universe {k = 1), we should make an analysis to see if the list of singularities given above is also valid in this case. 
It is straightforward to see that all the singularities listed above can be reproduced for a particular choice of the 
effective EoS. To show how the cosmic bounds behave for each type of singularity, we could write an explicit solution 
of the FRW equations, expressed as a function of time depending on free parameters that will be fixed for each kind 
of singularity. Then, the Hubble parameter may be written as follows 



16ttG I 
n(n — 1) ' 



g(*) = J „,„ ^ p-- s =H 1 (t a -t) m + H 0) (36) 



where m is a constant properly chosen for each type of singularity. Note that this is just a solution that ends in the 
singularities mentioned above, but there are other solutions which also reproduce such singularities. We will study 
how the cosmic bounds behave near each singularity listed above. 

As pointed out in Ref. @, [j| , around a singularity quantum effects could become important as the curvature of the 
Universe grows and diverges in some of the cases. In other words, approaching the finite-time future singularity the 
curvature grows and universe reminds the early universe where quantum gravity effects are dominant ones because of 
extreme conditions. Then one has to take into account the role of such quantum gravity effects which should define 
the behaviour of the universe just before the singularity. Moreover, they may act so that to prevent the singularity 
occurrence. In a sense, one sees the return of quantum gravity era. However, the consistent quantum gravity theory 
does not exists so far. Then, in order to estimate the influence of quantum effects to universe near to singularity 
one can use the effective action formulation. We will apply the effective action produced by conformal anomaly 
(equivalently, the effective fluid with pressure/energy-density corresponding to conformal anomaly ones) because of 
several reasons. It is known that at high energy region (large curvature) the conformal invariance is restored so one 
can neglect the masses. Moreover, one can use large N approximation to justify why large number of quantum fields 
may be considered as effective quantum gravity. Finally, in the account of quantum effects via conformal anomaly we 
keep explicitly the graviton (spin 2) contribution. The conformal anomaly Ta has the following well-known form 

Ta — b + Hqr^ + b'G + b"UR . (37) 

Here we assume for simplicity a 3+1 dimensional spacetime. Then, F is the square of a 4D Weyl tensor and G is the 
Gauss-Bonnet invariant, 

F=^R 2 - 2R ij R ij + Rij k iR m , G = R 2 - Mi, .If + R ijk iR ijkl . (38) 

The coefficients b and b 1 in (|3"7| are described by the number of N scalars, N x y 2 spinors, N\ vector fields, N2 gravitons 
and -/Vjjd higher derivative conformal scalars. They can be written as 

_ N + 6N l/2 + 12Ni + 6ILZV2 - 8iV HD , _ N + llN 1/2 + 62Ni + UIIN2 - 28N HD 
b = 120(4^)2 ' b = 360(4tt)2 ' (39) 

As b" is arbitrary it can be shifted by a finite renormalization of the local counterterm. The conformal anomaly Ta 
can be written as Ta = —Pa + Spa, where pA and pa are the energy and pressure densities respectively. By using 
([57)1 and the energy conservation equation pa + 3H(pa + Pa) — 0, one obtains the following expression for pa (5lll3|. 

Pa = — -7 I dta 4 HT A 



~ I dta A H 
a 4 



-!2bH 2 + 24b'(-H 2 + H 2 H + H 4 ) - (4b + 6b")(H + 7HH + 4H 2 + 12H 2 H) . (40) 
The quantum corrected FRW equation is given by 

H 2 = —(p + pa)- ^ ■ (41) 

We study now how the bounds behave around the singularity in the classical case when no quantum effects are added, 
and then include the conformal anomaly (1571) quantum effects in the FRW equations. We will see that for some cases 
the violation of the cosmic bound can be avoided. 
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A. Big Rip Singularity 



This type of singularity has been very well studied and has become very popular as it is a direct consequence in the 
majority of the cases when the effective EoS parameter is less than —1, the so-called phantom case [HQ- Observations 
currently indicate that the phantom barrier could have already been crossed or it will be crossed in the near future, 
so a lot of attention has been paid to this case. It can be characterized by the solution (|36|) with m < —1, and this 
yields the following dependence of the total energy density on the scale factor near the singularity, when a 3> 1, for a 
closed Universe (k = 1), 



P = 



n(n — 1) 
16ttG 



H z 



1(1+1 



for t — > t s 



(42) 



where we have chosen H\ = 2/n\l + w\ with w < — 1, m — — 1 and Hq = for clarity. This solution drives the 
Universe to a Big Rip singularity for t — >■ t s , where the scale factor diverges. If the singularity takes place, the bound 
(f3"2"j) has to be violated before this happens. This can be seen from equation (|4"2")l . as the Casimir energy behaves as 
Ec oc a"' 10 ' while the Bekenstcin-Hawking energy goes as Ebr oc a"~ 2 . Then, as w < —1, the Casimir energy grows 
faster than the BH energy, so close to the singularity where the scale factor becomes very big, the value of Eq will 
be much higher than Eb, thus violating the bound (1521 . Following the postulate from Ref.@ one could interpret the 
bound (|32l) as the limit where General Relativity and Quantum Field Theory converge, such that when the bound is 
saturated quantum gravity effects should become important. QG corrections could help to avoid the violation of the 
bound and may be the Big Rip singularity occurrence. As this is just a postulate based on the CV formula, which is 
only valid for special cases as shown in the sections above, the bound on Eq could not be valid for any kind of fluid. 

Let us now include the conformal anomaly (I37[) as a quantum effect that becomes important around the Big Rip. 
In such a case there is a phase transition and the Hubble evolution will be given by the solution of the FRW equation 
(I4ip . Let us approximate to get some qualitative results, assuming 3 + 1 dimensions. Around t s the curvature is large, 
and \pa\ >> {3/k 2 )H 2 + k/a 2 . Then p — -pA, and from (|lD])we get 



p + AHp = H 



\2bH 2 + 2ib'(-H 2 + H 2 H + iJ 4 ) - (46 + 66") (H + 7HH + AH 2 + \2H 2 H) 



We assume that the energy density, which diverges in the classical case, behaves now as 

P~ {t s ~t) x , 



(43) 



(44) 



where A is some negative number. By using the energy conservation equation p + 3iJ(l + w)p — 0, the Hubble 
parameter goes ss H ^ 1 / (t s — t) . We can check if this assumption is correct in the presence of quantum effects by 
inserting both results in Eq. (|43l) . We get 



P 



3H 4 (-m + 24b') 



(45) 



Hence as b > and b' < 0, p becomes negative, which is an unphysical result. Thus p should not go to infinity in the 
presence of the quantum correction. This is the same result as obtained in Ref. [5( where numerical analysis showed 
that the singularity is moderated by the conformal anomaly, so that the violation of the bound that naturally occurs 
in the classical case can be avoided/postponed when quantum effects are included. 



B. Sudden singularity 



This kind of singularity is also problematic with respect to the bounds, but as the energy density p does not diverge, 
the violation of the bound may be avoided for some special choices. The sudden singularity can be described by the 
solution (|36| with < m < 1, and constants Hq i > 0. Then the scale factor goes as 



a{i) oc exp 



m + 1 



(46) 



which gives a(t) ~ e ot (de Sitter) close to t s . From the first FRW equation the total energy density becomes 



H 2 (t) + 



[H^-tr+Ho] 2 +exp 



2H t 



(47) 
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which tends to a constant p ~ Hq + e~ 2H ° ts for t — > t s . Then the Casimir energy grows as Ec oc iJ^a™ + a"~ 2 , while 
Ebh oc a™ -2 close to t s . The BH energy grows slower than the Casimir energy and the bound is violated for a finite t. 
However, by an specific choice of the coefficients, the violation of the bound ([32]) could be avoided. For H = 0, and 
by some specific coefficients, the bound could be obeyed. In general, it is very possible that Eq exceeds its bound. In 
the presence of quantum corrections, the singularity can be avoided but the bound can still be violated, depending 
on the free parameters for each model. We may assume that in the presence of the conformal anomaly for n — 3, the 
energy density grows as [f| 

p = Po + P i(t s - t) x , (48) 

where p and p\ are constants, and A is now a positive number. Then the divergences on the higher derivatives of the 
Hubble parameter can be avoided, as is shown in Ref.Q. Nevertheless, Eq still grows faster than -Ebh, such that the 
Universe has to be smaller than a critical size in order to hold the bound (1321) as is pointed in Ref.[7[ for the case of 
a vacuum dominated universe. 



C. Type III singularity 



This type of singularity is very similar to the Big Rip, in spite of the scale factor a(t) being finite at the singularity. 
The solution (|36|) reproduces this singularity by taking — 1 < m < 0. The scale factor goes as 



a(t) = a s exp 



H 



1 



{ts - t) 



m+1 



where for simplicity we take Ho — 0. Then, for t —> t s , the scale factor a(t) 
singularity, let us write it in terms of the time instead of the scale factor, 

E c oc a s H 1 (t s - t) + a s , 



(49) 



To see how Ec behaves near the 



(50) 



where m < 0. Hence, the Casimir energy diverges at the singularity, while Ebh oc a"~ 2 takes a finite value for the 
singularity time t s , so the bound is clearly violated long before the singularity. Then, in order to maintain the validity 
of the bound ([32]). one might assume, as in the Big Rip case, that GR is not valid near or at the bound. Even if 
quantum effects are included, as was pointed in Ref.[5], for this type of singularity the energy density diverges more 
rapidly than in the classical case, so that the bound is also violated in the presence of quantum effects. 



D. Type IV singularity 



For this singularity, the Hubble rate behaves as 

H = H 1 {t) + {t s -t) a H 2 {t) . (51) 

Here H\ (t) and H2 (t) are regular function and do not vanish at t — t s . The constant a is not integer and larger than 
1. Then the scale factor behaves as 

lna(i)~ J dtHi[t) + J dt(t s - t) a H 2 (t) . (52) 

Near t — t s , the first term dominates and every quantities like p, p, and a etc. are finite and therefore the bound (|32[) 
would not be violated near the singularity. 



E. Big Bang singularity 



When the matter with w > coupled with gravity and dominates, the scale factor behaves as 

a^t^TTT^) , (53) 

Then there appears a singularity at t = 0, which may be a Big Bang singularity. Although the Big Bang singularity is 
not a future singularity, we may consider the bound (|32p when t ~ 0. Since n(l + w) > 2, the energy density behaves 
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as p ~ a~ n ( 1+w ^ and therefore the Casimir energy behaves as Eq ~ a - "™. On the other hand, we find £bh ~ a™ -2 . 
Then when n > 2 or when n > 2 and it) > 0, £"c dominates when a — > 0, that is, when t — > 0, and the bound (I3"2"j) is 
violated. This tells us, as expected, that quantum effects become important in the early universe. 

Above, we have thus explored what happens near the future cosmic singularities. We have seen that in general, 
and with some very special exceptions on the case of Type II and Type IV, the bound will be violated if one assumes 
the validity of GR close to the singularity. Even if quantum corrections are assumed, it seems that the bound will be 
violated, although in the Big Rip case the singularity may be avoided when quantum effects are incorporated. It is 
natural to suggest, in accordance with Verlinde, that the bound on the Casimir energy means a finite range for the 
validity of the classical theory. When this kind of theory becomes saturated, some other new quantum gravity effects 
have to be taken into account. We conclude that the universality of the bound ([521 is not clear and may hold just for 
some specific cases, like the radiation dominated Universe. 



IV. F(i?)-GRAVITY AND THE CARDY- VERLINDE FORMULA 



It is known that modified gravity (for general introduction, see [l4j ) may be presented in the form of generalized 
fluid with inhomogeneous EoS [Hj]. This type of theories which became popular recently may pretend to unify the 
early-time inflation theory with the theory describing the late-time acceleration (T^]. We specify here a modified F(R)- 
gravity modeled as an effective fluid and construct the corresponding CV formula for it. The action that describes 
i ? (i?)-gravity is given by 



S =^J d n+1 x^(F(R) + L m ) 



(54) 



where L m represents the matter Lagrangian and k 2 — 8irG. The field equations are obtained by varying the action 
(|54p with respect to the metric g^ v , 



R^F'(R) - \g„ v F{R) + g^DF'(R) - V M V„F'CR) = k 2 T^ 



(55) 



Here is the energy-momentum tensor for the matter filling the Universe, and we have assumed a 1 + 3 spacetime 
for simplicity. For closed 3 + 1 FRW Universe, the modified FRW equations are expressed as 

1 



-F(R) 



-F(i£) - 3(H 2 + H)F'(R) + 3HF"{R)R = n 2 p m , 
3H 2 + H+^ F'(R) - [(d tt F'(R)) + 2H(d t F'(R))} = K 2 p m , 



(56) 



where primes denote derivatives respect to R and dots with respect to t. These equations can be rewritten in order 
to be comparable with those of standard GR. For such a propose the geometric terms can be presented as an effective 
energy-density Pf(r) an d a pressure Pf(r), 



1 



1 



3F'{R)' 3F'(R) 



RF'(R) - F(R) 



3HRF"(R) 



2H + 3H H — - 



F'{R) 



F'(R) 



R 2 F'"{R) + 2HRF"(R) + RF"(R) + -(F(R) - RF'(R)) 



(57) 



Then, an EoS for the geometric terms can be defined as Pf(r) = w f(R)Pf(R) ■ We can define an effective energy-density 
P = Pm/F'(R) + Pf(r) and pressure p = p rn /F'(R) +Pf(r)- Hence, for some special cases the formula for the entropy 
developed in the second section can be obtained in F(i?)-gravity (for an early attempt deriving a CV formula in a 
specific version of i 7 '(i?)-gravity, see [I3])- For example, for an F(R) whose solution gives p oc a~ 3 ^ 1+Wc!! \ the formula 
for the entropy ([9]) is recovered although in general, as in the cases studied above, no such expression can be given. 
On the other hand, one could assume that the geometric terms do not contribute to the matter sector. Supposing 
a constant EoS matter fluid, the expression for the entropy is given by ©, although the cosmic Cardy formula (fTTj) 
has not the same form and in virtue of the modified first FRW equation ([56)) the form of the Hubble entropy Sh, 
the total energy E, and the Bekenstein energy -Ebh, will be very different. It is not easy to establish correspondence 
between two such approaches. Note that using the effective fluid representation the generalized CV formula may be 
constructed for any modified gravity. 



11 



Now we consider the case where F(R) behaves as 

F(R) ~ R a , (58) 
when the curvature is small or large. Then if the matter has the EoS parameter w > — 1, by solving (|56|) we find 



C t when , > 

a 2Q . (59) 

I (ts-t)— when^ y <0 

Then there may appear a singularity at t = 0, which corresponds to the Big Bang singularity, or at £ = t s , which 
corresponds to the Big Rip singularity. Since the Casimir energy behaves as Ec ~ a~ nw but -Ebh ~ a™~ 2 , only when 
£ — > 0, i?c dominates in case that n > 2 and w > or in case that rt > 2 and w > 0. Even in the phantom phase 
where n ^° w ^ < 0, the bound (|32j) is not violated. 

Let us now consider de Sitter space solution in F(i?)-gravity (for review of CV formula in dS or AdS spaces, see [18|). 
As was pointed in Ref. 19], almost every function F(R) admits a de Sitter solution. This can be easily seen from 
the first FRW equation in (f56|). A de Sitter solution is given by a constant Hubble parameter H(t) = Ho; then by 
inserting in (|56[) we obtain the following algebraic equation, 

Here i?o = 12i?o an d the contribution of matter is neglected. Then, for positive roots H of this equation, the 
corresponding F(R) leads to the de Sitter solution which may describe inflation or dark energy. In this case the 
formula for the entropy @ can be reproduced for w = — 1, and even the universal bound (|32p can hold by taking a 
critical size of the Universe. The formula that relates the cosmic bounds in (135[) is easily obtained also in F(i?)-gravity 
for a de Sitter solution. In such a case one can identify 

H V aV F(R ) V 

Sii -^G> 5B "24GFW BH "^' ( } 

which corresponds to the first FRW equation written as 5^ + (5b — 5bh) 2 = 5|. Thus, one can conclude that 
dynamical entropy bounds are not violated for modified gravity with de Sitter solutions. Note that quantum gravity 
effects may be presented also as an effective fluid contribution. In case when de Sitter space turns out to be the 
solution, even with the account of quantum gravity the above results indicate that dynamical cosmological/entropy 
bounds are valid. In other words, the argument indicates the universality of dynamical bounds. It seems that their 
violation is caused only by future singularities if they are not cured by quantum gravity effects. Note that a large 
number of modified gravity theories do not contain future singularities; they are cured by higher derivatives terms. 



V. DISCUSSIONS 



In summary, we have derived a generalized CV formula for multicomponent, interacting fluids, generalized in the 
sense that an inhomogeneous EoS (including viscous fluid) was assumed. We also considered modified -F(i?)-gravity, 
using its fluid representation. We showed that for some special cases the formula is reduced to the standard CV 
formula expressing the correspondence with 2d CFT theory. The dynamical entropy bound for all above cases was 
found. The universality of dynamical entropy bound near all four types of the future singularity, as well as the initial 
Big Bang singularity, was investigated. It was proved that except from some special cases of Type II and Type IV 
singularity the dynamical entropy bound is violated near the singularity. Taking into account quantum effects of 
conformally invariant matter does not improve the situation. 

One might think that the dynamical entropy bound is universal and that its violation simply indicates that the 
situation will be changed with the introduction of quantum gravity effects. However, arguments given below indicate 
that it is not the case and that the future singularity is the domain where all known physical laws and equations are not 
valid. Indeed, taking account of quantum effects such as done in section VI does not improve the situation with respect 
to non-universality of the dynamical entropy bound. From another side, it was shown that the dynamical entropy 
bound is valid for the de Sitter solution. Having in mind that quantum gravity corrections may always be presented 
as a generalized effective fluid, one sees that the dynamical cosmological bound is not valid near the singularity (even 
when account is taken of Quantum Gravity). It is only when modified gravity (with or without quantum corrections) 
is regular in the future, like the models of Refs. [3. [20L l2~j| where the future universe is asymptotically de Sitter, that 
the dynamical bound remains valid. Hence, the problem of non-universality of dynamical entropy bound is related to 
the more fundamental question about the real occurrence of a future singularity. It remains a challenge to find any 
observational indications for the structure of the future universe. 
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Appendix A: Account of Casimir effect in the CV formula 



A different way to account for the Casimir energy is to relate the effect to the single length parameter in the 
theory, that is the scale factor, and to assume that we can exploit the expression for the Casimir effect in a perfectly 
conducting spherical shell. The method consists essentially in identifying the scale factor with the radius of the shell. 
This proposal has been studied in Ref.|22J, and as shown, could have some important effects on the cosmological 
history. Here, we want to deduce the expression for the entropy following this alternative approach. The expression 
for the Casimir energy according to [22[ can be expressed as follows 

E ° = ffij • < A1 > 

where L is an auxiliary length that has been introduced due to the non-dimensional nature of a. This is the same 
form as the one encountered for a perfectly conducting shell (see Ref.[23j]), where it was found that 

C = 0.09235 . (A2) 

Hence, we assume that value of C is much less than the unity, which is physically reasonable in view of the conventional 
feebleness of the Casimir force. Also, for simplicity we shall assume a 1 + 3 FRW Universe. The expression (|A2p 
corresponds to a Casimir pressure 

-1 dE c 

Pc 8ir(La) 2 d{La) 87tL 4 a 4 ' 1 ' 

and leads consequently to a Casimir energy density p c cx 1/a 4 , which means that the Casimir fluid has an EoS 
parameter given by wq = 1/3 in order to obey the energy conservation equation Then, the Casimir energy 
density and pressure are given by 

Pc = ■ ( A4 ) 



^ c 87rL 4 a 4 ' KC 8nLW 

Now we assume that the Universe is filled with a perfect fluid p m = w m p rn , where p — E m /V. By following the 
rescaling properties ©, and assuming now w m = 1/3, the extensive and the Casimir energy are written as a function 
of the entropy, 

The total energy is E = E m + 1/2Eq- Then, the expression for the entropy is very similar to the one obtained in 
Ref.Q, 

2?tLa 



^Ec(2E-E c ) , (A6) 



where a and /3 are arbitrary constants. We see that this expression, except from the constants, is equal to (fTU)) ; the 
constants can be absorbed by af3. Then, this approach also supports the formula (|10p . but like that, it is not possible 
to extend the formalism to matter fluids with arbitrary w m . In that case the expression for the entropy would not be 
constant and the first law of thermodynamics would be violated. 

Hence, we have obtained a twofold desription of the Casimir effect. The Casimir energy appears explicitly in the 
entropy formula. Moreover, all quantities are constructed from an effective Casimir fluid. 



[1] R. R. Caldwell, Phys. Lett. B 545, 23 (2002) [arXiv:astro-ph/9908168]; 

R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 91, 071301 (2003) [arXiv:astro-ph/0302506]; 



13 



B. Mclnnes, JHEP 0208 (2002) 029 [arXiv:hep-th/0112066]; 

S. Nojiri and S. D. Odintsov, Phys. Lett. B 562, 147 (2003) [arXiv:hep-th/0303117]; 
S. Nojiri and S. D. Odintsov, Phys. Lett. B 571, 1 (2003) [arXiv:hep-th/0306212]; 
V. Faraoni, Int. J. Mod. Phys. D 11, 471 (2002) [arXiv:astro-ph/0110067]; 

P. F. Gonzalez-Diaz, Phys. Lett. B 586, 1 (2004) [arXiv:astro-ph/0312579]; TSPU Vestnik 44N7, 36 (2004) [arXiv:hep- 
th/0408225]; 

M. Sami and A. Toporensky, Mod. Phys. Lett. A 19, 1509 (2004) [arXiv:gr-qc/0312009]; 

H. Stefancic, Phys. Lett. B 586, 5 (2004) [arXiv:astro-ph/0310904]; 

L. P. Chimento and R. Lazkoz, Phys. Rev. Lett. 91, 211301 (2003) [arXiv:gr-qc/0307111]; 

Mod. Phys. Lett. A 19, 2479 (2004) [arXiv:gr-qc/0405020]; 

J. G. Hao and X. Z. Li, Phys. Lett. B 606, 7 (2005) [arXiv:astro-ph/0404154]; 

E. Babichev, V. Dokuchaev and Yu. Eroshenko, Class. Quant. Grav. 22, 143 (2005) [arXiv:astro-ph/0407190]; 
X. F. Zhang, H. Li, Y. S. Piao and X. M. Zhang, Mod. Phys. Lett. A 21, 231 (2006) [arXiv:astro-ph/0501652]; 

E. Elizalde, S. Nojiri, S. D. Odintsov and P. Wang, Phys. Rev. D 71, 103504 (2005) [arXiv:hep-th/0502082]; 
M. P. Dabrowski and T. Stachowiak, Annals Phys. 321, 771 (2006) [arXiv:hep-th/0411199]; 

F. S. N. Lobo, Phys. Rev. D 71, 084011 (2005) [arXiv:gr-qc/0502099]; 

R. G. Cai, H. S. Zhang and A. Wang, Commun. Theor. Phys. 44, 948 (2005) [arXiv:hep-th/0505186]; 

I. Y. Arefeva, A. S. Koshelev and S. Y. Vernov, Theor. Math. Phys. 148, 895 (2006) [arXiv:astro-ph/0412619]; Phys. Rev. 
D 72, 064017 (2005) [arXiv:astro-ph/0507067]; 

H. Q. Lu, Z. G. Huang and W. Fang, arXiv:hep-th/0504038; 

W. Godlowski and M. Szydlowski, Phys. Lett. B 623, 10 (2005) [arXiv:astro-ph/0507322]; 
J. Sola and H. Stefancic, Phys. Lett. B 624, 147 (2005) [arXiv:astro-ph/0505133]; 

B. Guberina, R. Horvat and H. Nikolic, Phys. Rev. D 72, 125011 (2005) [arXiv:astro-ph/0507666]; 

M. P. Dabrowski, C. Kiefer and B. Sandhofer, Phys. Rev. D 74, 044022 (2006) [arXiv:hep-th/0605229]; 
E. M. Barbaoza and N. A. Lemos, arXiv:gr-qc/0606084. 

D. Saez-Gomez, Gen. Rel. Grav. 41, 1527 (2009) [arXiv:0809.1311 [hep-th]]. 

E. Elizalde, S. Nojiri, S. D. Odintsov, D. Saez-Gomez and V. Faraoni, Phys. Rev. D 77, 106005 (2008) [arXiv:0803.1311 
[hep-th]]. 

[2] E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70, 043539 (2004) [arXiv:hep-th/0405034]. 
[3] J. D. Barrow, Class. Quant. Grav. 21, L79 (2004) [arXiv:gr-qc/0403084]; 

S. Cotsakis and I. Klaoudatou, J. Geom. Phys. 55, 306 (2005) [arXiv:gr-qc/0409022]; 

S. Nojiri and S. D. Odintsov, Phys. Lett. B 595, 1 (2004) [arXiv:hep-th/0405078]; S. Nojiri and S. D. Odintsov, Phys. 
Rev. D 70, 103522 (2004) [arXiv:hep-th/0408170]; 

J. D. Barrow and C. G. Tsagas, Class. Quant. Grav. 22, 1563 (2005) [arXiv:gr-qc/0411045]; 

M. P. Dabrowski, Phys. Rev. D 71, 103505 (2005) [arXiv:gr-qc/0410033]; A. Balcerzak and M. P. Dabrowski, Phys. Rev. 
D 73, 101301(R) (2006) [arXiv:hep-th/0604034] ; 

L. Fernandez- Jambrina and R. Lazkoz, Phys. Rev. D 70, 121503(R) (2004) [arXiv:gr-qc/0410124]; L. Fernandez- Jambrina 
and R. Lazkoz, Phys. Rev. D 74, 064030 (2006) [arXiv:gr-qc/0607073] ; arXiv:0805.2284 [gr-qc]; 

P. Tretyakov, A. Toporensky, Y. Shtanov and V. Sahni, Class. Quant. Grav. 23, 3259 (2006) [arXiv:gr-qc/0510104]; 

H. Stefancic, Phys. Rev. D 71, 084024 (2005) [arXiv:astro-ph/0411630]; 

A. V. Yurov, A. V. Astashenok and P. F. Gonzalez-Diaz, Grav. Cosmol. 14, 205 (2008) [arXiv:0705.4108 [astro-ph]]; 

I. Brevik and O. Gorbunova, Eur. Phys. J. C 56, 425 (2008) [arXiv:0806.1399 [gr-qc]]; 

M. Bouhmadi-Lopez, P. F. Gonzalez-Diaz and P. Martin-Moruno, Phys. Lett. B 659, 1 (2008) [arXiv:gr-qc/0612135]; 
arXiv:0707.2390 [gr-qc]; 

M. Sami, P. Singh and S. Tsujikawa, Phys. Rev. D 74, 043514 (2006) [arXiv:gr-qc/0605113]; 

C. Cattoen and M. Visser, Class. Quant. Grav. 22, 4913 (2005) [arXiv:gr-qc/0508045]; 
J. D. Barrow and S. Z. W. Lip, arXiv:0901.1626 [gr-qc]; 

T. Koivisto, Phys. Rev. D 77, 123513 (2008) [arXiv:0803.3399 [gr-qc]]; 
M. Bouhmadi-Lopez, Y. Tavakoli and P. V. Moniz, arXiv:091 1.1428 [gr-qc]. 
[4] S. Nojiri and S. D. Odintsov, Phys. Rev. D 78, 046006 (2008) [arXiv:0804.3519 [hep-th]]; K. Bamba, S. Nojiri and 
S. D. Odintsov, JCAP 0810, 045 (2008) [arXiv:0807.2575 [hep-th]]; 

S. Capozziello, M. De Laurentis, S. Nojiri and S. D. Odintsov, arXiv:0903.2753 [hep-th]. 

K. Bamba, S.D. Odintsov, L. Sebastiani and S. Zerbini, arXiv:0911.4390. 
[5] S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71 063004 (2005) [arxiv:hep-th/0501025] 
[6] E. Verlinde, arXiv:hep-th/0008140. 

[7] D. Youm, Phys. Lett. B 531, 276 (2002) [arXiv:hep-th/0201268]. 
[8] J. L. Cardy, Nucl. Phys. B 270 (1986) 186. 

[9] S. Nojiri and S. D. Odintsov, Phys. Rev. D 72, 023003 (2005) [arXiv:hep-th/0505215]; Phys. Lett. B 639, 144 (2006) 
[arXiv:hep-th/0606025] . 
[10] D. Saez-Gomez, Grav. Cosmol. 15 134 (2009), [arxiv:0804.4586 [hep-th]]. 

[11] I. Brevik, in Horizons in World Physics, Vol. 246 Quantum Cosmology Research Trends, edited by A. Reimer (Nova Science 

Publ., New York, 2005) p. 165 [arXiv:gr-qc/0404095] . 
[12] I. Brevik and S. D. Odintsov, Phys. Rev. D 65, 067302 (2002) [arXiv:gr-qc/0110105]; 

I. Brevik, Phys. Rev. D 65, 127302 (2002) [arXiv:gr-qc/0204021]; 



14 



I. Brevik, Int. J. Mod. Phys. A 18, 2145 (2003) [arXiv:gr-qc/0301026]. 
[13] S. Nojiri and S. D. Odintsov, Int. J. Mod. Phys. A 16, 3273 (2001) [arXiv:hep-th/0011115]. 
[14] S. Nojiri and S. D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007) [arXiv:hep-th/0601213]. 
[15] S. Capozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B 634, 93 (2006) [arXiv:hep-th/0512118]; 

S. Capozziello, S. Nojiri, S. D. Odintsov and A. Troisi, Phys. Lett. B 639, 135 (2006) [arXiv:astro-ph/0604431]. 
[16] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003) [arXiv:hep-th/0307288]. 

[17] I. Brevik, S. Nojiri, S. D. Odintsov and L. Vanzo, Phys. Rev. D 70, 043520 (2004) [arXiv:hep-th/0401073]. 
[18] R. G. Cai, Nucl. Phys. B 628, 375 (2002) [arXiv:hep-th/0112253]; Phys. Lett. B 525 331 (2002) [arxiv:hep-th/0111093]. 
[19] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov and S. Zerbini, JCAP 0502, 010 (2005) [arXiv:hep-th/0501096]. 
[20] M. C. B. Abdalla, S. Nojiri and S. D. Odintsov, Class. Quant. Grav. 22, L35 (2005) [arXiv:hep-th/0409177]. 
[21] A. Dev, D. Jain, S. Jhingan, S. Nojiri, M. Sami and I. Thongkool, Phys. Rev. D 78, 083515 (2008) [arXiv:0807.3445 
[hep-th]]; 

T. Kobayashi and K. I. Maeda, Phys. Rev. D 79, 024009 (2009) [arXiv:0810.5664 [astro-ph]]; 

I. Thongkool, M. Sami, R. Gannouji and S. Jhingan, Phys. Rev. D 80, 043523 (2009) [arXiv:0906.2460 [hep-th]]; 

S. Appleby, R. Battye and A. Starobinsky, arXiv:0909.1737 [astro-ph.CO]. 
[22] I. Brevik, O. Gorbunova and D. Saez-Gomez, arXiv:0908.2882 [gr-qc]; 

O. Gorbunova and D. Saez-Gomez, arXiv:0909.5113 [gr-qc]. 
[23] K. A. Milton, L. L. DeRaad and J. Schwinger, Annals Phys. 115, 388 (1978). 



